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Abstract 

In this work, we analytically investigate a degenerating PDE system for phase separation 
and complete damage processes considered on a nonsmooth time-dependent domain with mixed 
boundary conditions. The evolution of the system is described by a degenerating Cahn-Hilliard 
equation for the concentration, a doubly nonlinear differential inclusion for the damage variable 
and a degenerating quasi-static balance equation for the displacement field. All these equations 
are highly nonlinearly coupled. Because of the doubly degenerating character of the system, the 
doubly nonlinear differential inclusion and the nonsmooth domain, the structure of the model is 
very complex from an analytical point of view. 

A novel approach is introduced for proving existence of weak solutions for such degenerating 
coupled system. To this end, we first establish a suitable notion of weak solutions, which consists 
of weak formulations of the diffusion and the momentum balance equation, a variational inequal¬ 
ity for the damage process and a total energy inequality. To show existence of weak solutions, 
several new ideas come into play. Various results on shrinking sets and its corresponding local 
Sobolev spaces are used. It turns out that, for instance, on open sets which shrink in time a 
quite satisfying analysis in Sobolev spaces is possible. The presented analysis can handle highly 
nonsmooth regions where complete damage takes place. To mention only one difficulty, infinitely 
many completely damaged regions which are not connected with the Dirichlet boundary may 
occur in arbitrary small time intervals. 

Key Words: Cahn-Hilliard system, phase separation, complete damage, elliptic-parabolic degen¬ 
erating systems, linear elasticity, weak solution, doubly nonlinear differential inclusions, existence 
results, rate-dependent systems. 

AMS Subject Classifications: 35K85, 35K55, 49J40, 49S05, 35J50, 74A45, 74G25, 34A12, 82B26, 
82C26, 35K92, 35K65, 35K35; 

1 Problem description 

Phase separation and damage processes occur in many fields, including material sciences, biology and 
chemical reactions. In particular, for the manufacturing and lifetime prediction of micro-electronic 
devices it is of great importance to understand the mechanisms and the interplay between phase- 
separation and damage processes in solder alloys. As soon as elastic alloys are quenched sufficiently, 
spinodal decomposition leads to a fine-grained structure of different chemical mixtures on a short 
time-scale (see |DM01j for numerical simulations and experimental observations). The long-term 
evolution is determined by a chemical diffusion process which tends to minimize the bulk and the 
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surface energy of the chemical substances. J.W. Cahn and J.E. Hilliard developed a phenomenolog¬ 
ical model for the kinetics of phase-separation in a thermodynamically consistent framework known 
as the Cahn-Hilliard equation [CH58| . for which an extensive mathematical literature exists. An 
overview of modeling and analytical aspects of the Cahn-Hilliard equation can be found in [E1189| . 
The recent literature is mainly focused on coupled systems. Eor instance, physical observations 
and numerical simulations reveal that mechanical stresses influence the developing of shapes of the 
chemical phases. A coupling between Cahn-Hilliard systems and elastic deformations have been 
analytically studied in [GarOni lBCD'*~n2l ICMPODl IGarn5al IGarOhbl IBPn51 IPZn8| . For numerical 
results and simulations we refer to |Wei01l IMerOSl IBB991 IGRWOll IBMIO] . Phase separation of the 
chemical components may also lead to critical stresses at phase boundaries due to swelling which 
result in cracks and formation of voids and are of particular interest to understand the aging process 
in solder materials, cf. [HGW9HIUSG071 IGUE+OTI IFKn9j . A fully coupled system consisting of the 
Cahn-Hilliard equation, an elliptic equation for the displacement field and a differential inclusion 
for the damage variable has been recently investigated in |HKllllHK13| . However, in |HKllllHK13| 
and in the most mathematical damage literature |BS04[ IGiaOSl IMR06LIMTIOI IKRZllj , it is usu¬ 
ally assumed that damage cannot completely disintegrate the material (i.e. incomplete damage). 
Dropping this assumption gives rise to many mathematical challenges. Therefore, global-in-time 
existence results for complete damage models are rare. Modeling and existence of weak solutions 
for purely mechanical complete damage systems with quasi-static force balances are studied in 
|BMR09l IMielll IHK12] and with visco-elasticity in [MRZlOl IRR12] . 

The main goal in the present work is to prove existence of weak solutions of a system coupling 
a differential inclusion describing damage processes with an elastic Cahn-Hilliard system in a small 
strain setting as in |HK11| but allowing for complete damage and degenerating mobilities with 
respect to the damage variable (see Theorem 14.5p . Note that the mobility is still assumed to be 
independent of the concentration variable (see |EG96j for an analysis of Cahn-Hilliard equations 
with degenerating and concentration dependent mobilities). The elasticity is considered to be linear 
and the system is assumed to be in quasi-static mechanical equilibrium since diffusion processes take 
place on a much slower time scale. The main modeling idea in [HK12] has been to formulate such 
degenerating system on a time-dependent domain which consists of the not completely damaged 
regions that are connected to the Dirichlet boundary. In this context, the concept of maximal 
admissible subsets is introduced to specify the domain of interest. 

In the following, we fix a bounded C^-domain G C M” and a Dirichlet boundary part D C cIG 
with > 0. 

A relatively open subset G of G is called admissible with respect to an '^"“^-measurable part 
D of the boundary 9G if every path-connected component Pq of G satisfies n D) > 0, 

where denotes the (n — l)-dimensional Hausdorff measure. The maximal admissible subset 

of G is denoted by 21^(0). With the notion of maximal admissible subsets, we can formulate our 
evolutionary system with a time-dependent domain. 

Degenerating PDE system on a time-dependent domain. 

Find functions 

ceC‘^{F]R), u€CI{F-,R^), zeC'^{Th;R), /ieC^(F;M) 

with the time-dependent domain F = {{x,t) G Gy | x G F{t)} and F{t) = 2t£)({z(t) > 0}) such that 


2 


























































the PDE system 


0 = diY {W^e{c,e{u),z)), 
ct = div(m(2;)V^), 

// = -Ac + ^,c(c) + W^cic, e{u), z), 

zt + C- div(|V2;|^“^Vz) + W^z{c, e{u),z) + f{z) = 0, 


is satisfied pointwise inmt(F) with the initial-boundary conditions 


c(0) = c°, 2:(0) = z^ 
u{t) = b{t) 

W^eic{t),e{u{t)),z{t)) • 1/ = 0 

z{t) = 0 

Vz{t) ■ ly = 0 

Vc{t) - 1 ^ = 0 

m{z{t))S/fi{t) ■ u = 0 


in F{0), 

on ri(t) := F{t) n D, 
on r 2 {t) := F{t) n {dn\D), 
on Fsit) :=dF{t)\F{t), 
on ri(t) U r2(i), 

on ri{t) u r2(t), 

on Ti{t) U r2(t). 


Here, C^{F; denotes the space of two times continuously differentiable functions with respect 
to the spatial variable x on the set F. 

The solution of the PDE system can physically be interpreted as follows: c denotes the chemical 
concentration difference of a two-component alloy, e{u) := ^(Vu+(Vm)^) the linearized strain tensor 
of the displacement u, z the damage profile describing the degree of damage (i.e. z = 1 undamaged 
and z = 0 completely damaged material point) and /r the chemical potential. Moreover, W denotes 
the elastic energy density, 'h the chemical energy density, / a damage dependent potential, m the 
mobility depending on the damage variable z and b the time-dependent Dirichlet boundary data 
for D. Capillarity effects are modeled by concentration gradients in the free energy which provokes 
the chemical potential to depend on the Laplacian Ac (see |Gur89j h 

We assume the following product structure for the elastic energy density: 

W{c,e,z) = g{z)ip{c,e) (1) 

with a non-negative monotonically increasing function g € C^([0,1]; R'*') such that the complete 
damage condition g(0) = 0 is fulfilled. The second function G C^(R x R”y^;R+) should have the 
following polynomial form 

(f{c, e) = ip^e : e + : e + ip^{c) (2) 

for coefficients G T(R”y^) with ip^ positive definite, ip"^ G C^(R;R”y^) and (p^ G C^(R). Note 
that homogeneous elastic energy densities of the type 

W (c, e, z) = ^zC(e - e*(c)) : (e - e*(c}) 

with eigenstrain e* depending on the concentration and stiffness tensor C are covered in this ap¬ 
proach. We remark that this approach allows to incorporate swelling phenomena of the chemical 
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species which, in turn, can lead to critical stresses for the initiation of damage. The mobility 
m € C([0,1];]R“*’), on the other hand, should satisfy the following condition for degeneracy, i.e. 

m{z) = 0 if and only if z = 0. (3) 


Plan of the paper 

A weak formulation of the degenerating system is given in Section [2] while the proof of existence 
is carried out in Section [3] and Section HI In the first step of the proof, a degenerated limit of the 
corresponding incomplete damage system coupled with elastic Cahn-Hilliard equations is performed 
(see Section ED- Due to the additional coupling the passage to the limit becomes quite involved and 
a conical Poincare inequality is used to control the chemical potential in a local sense. 

In the next step, we take material exclusions into account (see Section H]). They occur when not 
completely damaged components become disconnected to the Dirichlet boundary. The exclusions 
lead to jumps in the energy inequality. By a concatenation property and by Zorn’s lemma, we are 
able to prove the main theorem in this paper: a global-in-time existence result (see Theorem 14.51) . 
In addition, several T-results and various properties on shrinking sets and their corresponding local 
Sobolev spaces are used. 

2 Assumptions and notion of weak solutions 

As before, D C M” denotes a bounded C^-domain and D C dO, the Dirichlet boundary with 
> 0. Furthermore, T > 0 denotes the maximal time of interest and p > n a, fixed 

exponent. 

The weak formulation of the PDF system presented in this work will be based on an energetic 
approach and uses the associated free energy. Let G C D be a relatively open subset. Then, the 
free energy on G is given by 

£G{c,e,z) := -F T(c) + W{c,e,z) + f{z)j dx 

for c G H^G), e G L2(G;M”/^") and ^ G W^’P{G). 

We suppose the structural assumptions ([T|) and ([2]). Here, T G G^(R) denotes the chemical 
energy density, / G G^([0,1]) a damage dependent potential and g G G^([0,with gr(0) = 0 a 
function modeling the influence of the damage on the material stiffness. We assume the following 
growth conditions: 


\p\c)\,\ip%ic)\<C{l + \c\), (4a) 

\P^ic)\,\p%ic)\<C{l + \c\‘^), (4b) 

|'k,c(c)|<G(l + |c|2V2), (4c) 

r] < g'{z). (4d) 


The constants r],G > 0 are independent of c and z, and 2* denotes the Sobolev critical exponent. 
In the case n = 2, has to satisfy an r-growth condition for a fixed arbitrary r > 0 whereas we 
have no restrictions on in tbe one-dimensional case. 

The functions m G C([0,1];M^) will denote the mobility of the Cahn-Hilliard diffusion process 
and is assumed to satisfy ([3|). 
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Before we state the precise notion of weak solutions, we introduce some notation. 

The set {/ > 0} for a given function / £ has to be read as {x £ O | f{x) > 0} by 

employing the embedding ^ C(n) (because of p > n). The symbol Hg{G) denotes the 

subspace of H^{G) with zero trace on S', where G C M” is open and S is an 7^”“^-measurable subset 
of dil. 

Let / C M be an open interval. The subspace SBV{I] X) C BV[I]X) of special functions of 
bounded variation is defined as the space of functions / £ BV{I]X) where the decomposition 




for an f £ L^{I;X) exists. The function f is called the absolutely continuous part of the differential 
measure and we also write dff and Jf denotes the jump set. If, additionally, dff £ LP(/; X), p > 1, 
we write / £ SBVP{I;X). 

We say that a relatively open subset F C n-r is shrinking iff F{t) C F{s) for all 0 < s < t < T, 
where F{t) denotes the t-cut of F. For the analysis of our proposed system, we introduce local 
Sobolev functions. The space-time local Sobolev space R'^) for a shrinking set F is 

given by 


■■={v.F^R^ 


Vt£ (0,r], VC/CCF(t)open: u|c/x(o,t) € ^^(O, t; R^^))}. 


In the following, a weak formulation of the system above combining the ideas in [HKll] and 
[HK12| is given. 


Definition 2.1 (Weak solution of the coupled PDE system) A quadruple {c,u,z,p,) is called 
a weak solution with the initial-boundary data {c^,z^,b) if 

(i) Spaces: 

c e L°°{0,T; H\n)) n {H\n)y), c(0) = c°, 

u £ LfHl i^^{F; R”), u = b on Dt n F, 

z G L’=^{0,T-,W^’P{n)) n SBV\0,T-,L'^{n)), z{0)=z^, z+(t) = z-(t)l^p) 

p £ 

with e := e{u) £ L^(F;R"y^) where F := 21d({2:“ > 0}) F is a shrinking set 
(z'^{t) := Iims 4 ,i 2 ;(s) and z~{t) := lim^it 2 ;(s) in L?‘{Q)). 

(a) Quasi-static mechanical equilibrium: 

0=/ We(c(t),e(t),2;(t)) : e(C)dx (5) 

dm 

for a.e. t G (0,r) and for all ( G 
(Hi) Diffusion: 

/ dtC{c —c^)dxdt = / m{z)Xfj, ■ V C dx dt ( 6 ) 

J Q.'j' F 
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for all C € L‘^{0,T;H^{n)) n H\0,T-, L‘^{n)) with C{T) = 0 and 


j /iCdxdt = J (Vc-Vf +'i>^cic)C+ W^cic,e,z)(^ dxdt 


( 7 ) 


for all ( € L‘^{0,T-, H^{n)) with supp(C) C F. 

(iv) Damage variational inequality: 

0 < [ (\Vz{t)r^Vz{t) • VC+{W,Mt),e{t),z{t))+fiz{t)) + dfz 
JFit) ^ 

0 < z{f), 

0 > dfz{t) 

for a.e. t G (0,T) and for all f G with C < 0. 

(v) Energy inequality: 

£{t) + J{0,t) + 7>(0,t) < ^+(0) + Wext(0,t), 
where the terms in the energy inequality are as follows: 


dx (8) 


(9) 


• Energy: 

• Energy jump term: 


• Dissipated energy: 

• External work: 


£{t) := £F(^t){c{t),e{t),z{t)), 

j{o,t)-.= {£-is)-£^is)), 

s&Jzn{o,t] 

£~{t):= lim f ess inf 

s^t- V rG(s,t) / 

£+{t) := e+ 

with 0 < et < inf £F(t)ic{t),e{b{t) + f), z{t)), 


V{0,t) 



dx dt, 


Wext(0,t) := 


0 JF(s) 


W^eic,e{u), z) : €{dtb) dx dt. 


Remark 2.2 Under additional regularity assumptions, a weak solution satisfies 
lim ^^(^)(c(r),e(r),z(r)) = lim ess inf £'j7’(^)(c(??), e(??), 2 :(??)), 

T^S- T^S- 1?G(t,s) 

^l™^^^F(r)(c(r),e(r),z(r)) = c+ 

and reduces to the pointwise classical formulation presented in Section\^ (cf. IHK121 Theorem 3.7]). 

As mentioned before, the main aim of this paper is to prove global-in-time existence of weak 
solutions. Due to the analytical quite involved structure, we introduce the notion of weak solutions 
with a given fineness constant ry > 0. 
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Definition 2.3 (Weak solution with fineness rj) A tuple (c, e, u, z, p) together with a shrinking 
set F C Qrp ig called a weak solution with fineness rj > 0 if 

ce c(0) = c°, 

u £ u = b onDTH^niF), 

z £ L°°{0, T- W^’P{n)) n SBV^iO, T- L^{VL)), z ( 0 ) = z+{t) = z-{t)lF(t) 

// e 

such that e = e(it) in ^d{F) and the properties (ii)-(v) of Definition \2.1\ and 

VtG [0,r] : 2tz)({z-(t) >0}) CF{t) and C^{F{t)\^D{{z-{t) >0}))) < ry, 

yt£[0,T]\ U + ^Di{z-it)>0}) = F{t). 

t€C^* 

are satisfied. 

Here, Cz* denotes the set of cluster points from the right of the jump set Jz* of the function 
z* £ T; L^(H)) given by z*{t) := ) i-®-) is the damage profile of z 

restricted to the time-dependent domain ^£)({z~ > 0}) (1a : X —)• {0,1} denotes the characteristic 
function of a set A Q X). 

Remark 2.4 (i) If a weak solution {c, e,u, z, p,) on F with fineness rj according to Defi,nition \2.S\ 

satisfies Cz* = 0 then {c,u,z,p) is a weak solution according to Definition \2.1\ . 

(ii) If the initial value for the damage profile contains no completely damaged parts on H, we are 
also able to establish maximal local-in-time existence of weak solutions. More precisely, there 
exist a maximal value T > 0 with T < T and functions c, u, z and p defined on the time 
interval [0,T] such that {c,u,z,p) is a weak solution according to Definition \2. A Therefore, 
ifT <T, {c,u,z,p) cannot be extended to a weak solution on [0,T -|-e]. 

The remaining part of this paper is devoted to establish a global existence result in the sense of 
Defi n ition 12.31 fsee Theorem 14.5p . 

3 Degenerate limit of the regularized system 

In this section, we will start with a corresponding incomplete damage model coupled to an elastic 
Cahn-Hilliard system and then perform a limit procedure. For each e > 0, we define the regularized 
free energy as 

f,(c, e, z) := ^ [^\Vz\P + ^iVcp + T(c) + iy^(c, e, z) + /(z)) dx 

for functions c £ II^(fl}, e £ L^(D;M”y^), 2 £ The regularized elastic energy density and 

mobility are given by 

W^c, e, z) := {g{z) + e)(^(c, e), 
nf{z) := m{z) + e. 

A modification of the proof of Theorem 4.6 in |HK11] yields the following result. 
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Theorem 3.1 (e-regularized coupled PDE problem) Let e > 0. For given initial-boundary 
data Cg G H^{Q), G and G T; M"')) there exists a quadruple = 

{cs,U£, Zs, such that 

(i) Spaces: 

Ce e L°-{o,T-,H\n))nH\o,T;{H\n)r), cM = cl 

Ue G L°°(0,r; M")), = be on Dt, 

ZeeL^{0,T-,W^’P{n))nH\0,T-,L\Q)), z,(0) = 

yieeL\0,T-,H\n)) 


( 11 ) 


dx (12) 


(a) Quasi-static mechanical equilibrium: 

[ W^eic£{t),e{ueit)),Ze{t)) : e{C)dx = 0 (10) 

Jn 

for a.e. t G (0,T) and for all f G 
(in) Diffusion: 

/ (ce — c°)dtCdx dt = / m^(z£)V/i£ • V^dx dt 

J Sly J Q^t 

for all C G L^(0,T;//i(0)) with dtC G L‘^{Qt) o-nd C{T) = 0 and 

[ /^£(i)Cdx= f (\/Ce{t) -VC+ d>^e{c£{t))C+ W),{Ce{t),e{Ue{t)),Z, 

Jn Jn ^ 

for a.e. t G (0,T) and for all C G 

(iv) Damage variational inequality: 

0 < ^ (^\VZeit)f~‘^VZeit) ■ VC + {W(^{Ce{t),€{Ueit)),Zeit)) + f{Ze{f)) + dtZe{t) + r, 

0 < Zeif), 

0 > dtZe{t) 

for a.e. t G (0,T) and for all C G with C < 0, where re G L^{Dt) is given by 

re = -Xe{W(^{Ce,e{Ue),Ze) d- fiz))^ (14) 

with Xe G L°°{D) fulfilling Xe = d on {ze > 0} and 0 < Xe < 1 on {ze = 0}. 

(v) Energy inequality: 

£e{Ce{t),e{Ue{t)),Ze{t)) + (m%Ze)\VHel"^ + {diZel"^^ dx ds 

< Ee{Ce^e{Ue)^^e) + [ W^eiCe, e{Ue), Ze) : e{dtbe) dx ds (15) 

J nt 

holds for a.e. t G (0,T) where Ue minimizes 8e{c^, ef), z^) in with Dirichlet data 

6° := be(0) on D. 


dx 

(13) 
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Proof. The existence theorem presented in [HKllj can be adapted to our situation by considering 
the viscous semi-implicit time-discretized system (in a classical notation; we omit the e-dependence 
in the notation for the discrete solution at the moment): 


0 = div(W|(c^e(u^),z'=)) +<5div(|Vu'=|2Vu'=), 

^ = div(m"(/-i)V/), 

T 

= -Ac’^ + T c(c^) + + 6^—^— 

’ ’ T 

+ -e + C = div(|Vz^|P-2Vz*^) + e(u*^), /) + f{z^), 

'T ’ 


with the sub-gradients ^ € dI(^_oofi]{{z^ — z^~^)/t), C, S '9/[o^oo)(' 2^) and the discretization fineness 
r = T/M for M G N. The discrete equations can be obtained recursively starting from {c^,u^,z^) 
with vP := argmin„gj|^i(-f^.]un)^„l^=ftO|^£’£(c'^, u, by considering the Euler-Lagrange equations of 
the functional 


K^(c,u, z) := £s{c,u, z) + / -|VnEdx 

Jn 4 

,A :-1 2 


T 

+ 2 


^ — z 


T 


+ 


L^{n) 


c — c 


.k-l 


+ 6 




c — c 


.k-l 


L^{n)^ 


dehned on the subspace of H^{Q,) x lT^’^(n;M"') x with the conditions u\d = b{kT)\j:,, 

f^(c — c°) dx = 0 and 0 < z < z^~^ a.e. in n. The weighted scalar product (•, ■)x{z^-'^) is given by 


{u,v)x(,k-iy.= (^m%z^ ^)VA ^u,VA 

with the operator ^4 : Vq ^ Au := {nf{z^~^)SIu,V-) and the spaces 

Eo := = 

Vo-.= {ce{H\n)r\ (C,i)(hi)*xhi=o}. 

After passing the discretization hneness to 0, i.e. r —?> 0+, we obtain the corresponding regular¬ 
ized equations and inequalities for (fT0]l - (fT^ . A further passage 6^0'^ yields a weak solution as 
required. □ 


In the following, we always consider a given countable sequence —)■ 0^ as A; —>■ oo and omit 

the subscript k by simply writing e —>■ 0^. 

The existence proof is based on T-limit techniques of the reduced energy functionals of and 
with p£{c, e, z) := iy^(c, e, z) dx in order to gain a suitable energy estimate in the limit e —>■ 0+. 
The reduced energy functionals Cie and Jg are defined by 

{ min ^^(c, e(^ + C), z) if 0 < z < 1, 

CG//l,(n;K") 

oo else, 
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{ min -7>(c, e(^ + C), - 2 ) if 0 < 2 < 1, 

00 else. 

The T-limits of Cie and 5'e as e —5- O’*" exist in the topological space x lT^’°°(n; M"') x lTw^(17) 

and are denoted by £ and respectively. Here, denotes the space with its weak 

topology. The limit functional ^ is needed as an auxiliary construction in the following because 
it already captures the essential properties of <B. In the next section, we are going to prove some 
properties of the T-limit G: which are used in the global-in-time existence proof. 

Let (Cg, 6°, z^) —>■ (c°, b^, z^) as e —>■ 0^ be a recovery sequence for Gie A- G:. In particular, c® ^ 
in —>■ b^ in IT^’°°(H;M”) and ^ z^ in IT^’^(H). Furthermore, we set b^ := b — b^ + b^. 

For each e > 0, we obtain a weak solution {cs,U£, z^, He) for {c^,z'^,bs) according to Theorem 13.11 
By means of the energy estimate m, we deduce the following a-priori estimates for our proposed 
model: 

Lemma 3.2 There exists a constant C > 0 independent of £ such that 
(i) suptg[o,T] l|c£(t)||Hi(o) < C, 

l|6e||L2(Q,j,.]Knxn) < C with Cg := e£l{2^>o}: 

(Hi) suptg[o,T] lk£(t)||w^i,p(o) < C, 

(iv) \\^tZ£\\L2^^^^) < C, 

(v) \\W^{Ce,ee, Ze)\\L^{Q^T-,L^{n)) < C with Ce := e{Ue), 

(vi) \\mf {Ze)^/^VHe\\mnr-m ^ 

(vii) ||9tC£||i2(0,r;(Hi(o))*) < \\'m%Ze)VHe\\mnT-,R^) < C. 

Proof. Applying Gronwall’s lemma to the energy estimate (|15p and noticing the boundedness of 
£e:{c^, e{u^), z')) with respect to e G (0,1) show (iv), (vi) and 

Se{Ce{t),eeit),Ze(t)) < C (16) 

for a.e. t G (0,r) and all e G (0,1) and in particular (v). Item (vii) follows from (vi) and the 
diffusion equation (jlip . Taking the constraint J^Ce{t)dx = const and the restriction 0 < < 1 

into account, property (unp gives rise to 


l|c£||L°°(0,T;/fl(*^)) — 

lk£||L°°(0,T;VKl’P(O)) ^ C. 

Together with the control of the time-derivatives (iv) and (vii), we even obtain boundedness of 

l|c£(i)||Hi(o) < C, 
lke(^)llvKi>P(0) < C 

for every t G [0,T] and e G (0,1). Hence, (i) and (iii) are proven. 
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It remains to show (ii). By using the estimates 


0 < y. < 1 

0 < WUc^,e^,Zs) 


a.e. in Qx, 
a.e. in Qx, 


and the definition of in (|14p . we obtain (here (•)^ := max{-,0} and (•) := mm{-,0}) 

/ {W^ziCs,ee,Ce) + f'ize)) dxdt 

J {ze>0} 

= / {W^ziCe,es,Ce) + f'iZe)) dxdt 

J 

-[ {W^z{ce,es,Ci;) + f'{zs))^ dxdt - f (iyj.(c£, e^, Cg) +/'(z^))” dx dt 

J {ze=0} J {ze=0} 

< / {W^z{ce,es,Cs) + f'izs)) dxdt 

J 

-[ Xe{W^ziCe,ee,Ce) + f{Ze))^ dxdt- [ {f{Ze))~dxdt 
J {z^=0} J {Zi;=0} 

= [ {W^z{Ce,es,Ce)+fiZe)+re)dxdt- [ {f{Ze))~ dxdt. 

JQ't Ju^=o} 


IQt J{ze=0} 

To proceed, we test inequality (fT3]l with (^ = — 1 and integrate from t = 0 to t = T: 

[ (^, 2 (C£, ee,Ze) + f'{ze) + r^) dx dt < - I dtZe dx df. 

J J 

Combining ()17p and (dD, we end up with 

/ {W^z(ce,es,Ce) + f'izi;)) dxdt < - [ dtz^dxdt- f (/'(zg))" dxdf. 
J^Zr>0\ J Jiz^=0\ 


(17) 


(18) 


The right hand side is bounded because of the boundedness properties (hi) and (iv). But we also 
have (C > 0 and C > 0 are constants) 

/ {W^z{Ce,ee,Ce) + f'{Ze)) dxdt> [ (c\ee\‘^-C\Cs\‘^ + f'{Ze))dxdt 

J{ze> 0 } J{z^> 0 } ^ 

> C \ lei-P dxdt 


f leepdxdtT / ( - Clcep +/'(Ze)) dxdh 

Jut JU^>0} ^ ^ 


I^ZLt j {ze>0} 

Since the left hand side is bounded and dx dt is also bounded, we obtain 

(ii). ' □ 

We proceed by deriving convergence properties for c^, he, as e —>■ 0 . 

Lemma 3.3 There exists functions 

c€ L’=^{0,T-,H\n))nH\0,T;{H\n))*) with c(0) = c°, 
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z <0, dtz < 0 


e G 

2 € L°°{0,T;W^’P{n)) n H\0,T-, L^{n)) with ^(0) = 2 °, 
and a subsequence (we omit the index) such that for e —>■ O"'' 

(a) Ce^cin H\0,T-, {H^Q))*), 

Ce ^ c in T ; L^{Q,)) for all k > 1 , 1 < r < 2*, 

Cs{t) c{t) in H^{Q) for all t, 

Ce ^ c a.e. in Qt, 

(b) Ze^z in Fi(0,r;L2(f^)), 

Ze^ z in L'^(0,r; VFi’P(Q)) for all k>1, 

Ze{t) z{t) in W^’P{VL) for all t, 

Ze ^ z in C(Ot); 

(c) be ^ bin 

(d) Ce^ein 

Wfe{ce,ee,Ze) ^ We(c,e,z) in L^{{z > 

W%{ce,ee,Ze) ^ 0 in L‘^{{z = 

Wl{ce,ee,Ze) ^ W^c^e^z) in L^iz > 0};M"X-), 

Wfics,ee,Ze) ^0 in L^{{z = 0};]R->'-). 

Proof. The a-priori estimates from Lemma 13.21 and classical compactness theorems as well as 
compactness theorems from Lions and Aubin |Sim86| yield the convergence properties of (a) and 


★ 

Ze Z 

in L°°(0,r;LF^’P(Li)), 

(19a) 

Ze^ Z 

in H\0,T-,L^{n)), 

(19b) 

Ze^ Z 

in C{Qt)i 

(19c) 

Cg ^ e 

in L2(Llr;K^''''), 

(19d) 

lT_g(Cg, Be, Ze) ^ We 

in L2(Llr;K^'''') 

(19e) 


as e —)■ O'*' for a subsequence (we omit the subscript) and appropriate functions We, e and z. 

The proof of the strong convergence of Vzg in essentially relies on the uniform 

p-monotonicity property 


Cuc\x - y\P < (|xP ‘^x - \y\P ^y) ■ {x - y) (20) 

and on an approximation scheme {Cg} C LP(0, T; lT^’P(f2)) with Ce ^ 0 and 

Ce ^ 2 in LP(0,r;iyi’P(O)) as e ^ 0+, (21a) 

0 < Ce < a.e. in PIt for all e G (0,1). (21b) 

For details of the construction of (e, we refer to |HKllj . By using (|20p . we obtain the estimate: 
Cue / \Vze-Vz\Pdxdt< j {\Vze\^~‘^Vze-\Vz\P-‘^Vz)-V{ze-z)dxdt 
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\Vz,\P-^Vz,-V{ze-Q)dxdt 


+ / \VZe\P~‘^VZe-V{Ce-z)-\Vz\P-^Vz-V{z,-z)dxdt. 

J six’ 


B, 


The variational inequality (fT^t) tested with C(t) = Ce{t) — Zeit) and integrated from t = 0 to t = T 
yield 


Ae < 


< 


e{Ue), Ze) + f {Ze) + dtZ^ (Ce - Ze) dx dt 

/ 9'{Ze)(p^€{Ue) : €(Us){Ce - Zs)dxdt+ (Cs) : €{Ue){Cs - Ze) dx dt 

J fix' 

+ [ {Ce) + f' (Ze) +dtZe) iCe - Ze)dxdt. 

J Qx 


( 22 ) 


Here, we have used the fact that = 0 a.e. in {ze > 0} and Q — Ze = d a.e. in {ze = 0} by (I21bh . 
We treat each term on the right hand side of (1221) as follows: 


Due to (|4dl) and (I21bp . we obtain 

g'{ze)^^e{ue) : e(n£)(C£ - < 0 a.e. in Dt- 

Therefore, the hrst integral on the right hand side of (|22p is less or equal 0. 

Due to the estimates (^al) and (I21bp . and the estimates 

< Cg{ze), 

(follows from (I4dp and 5^(0) > 0) and 

9 iZe)\ee\^ < C{W‘^{Ce,ee,Ze) + \Cef) 

(follows from ([2]), ()4ap and ()4ap b the second integral on the right hand side of 
estimated as follows: 


can be 


/ 

J Qx 


V9^(Ce) : e(Ue)(Ce - Ze) dx dt 


^ ||v^ (Ce)||^oo(0^T;I,2(O;R»ixn))||e(^e)'\/iCe ^IL°°(0,T;I,2(Q;Rnxn)) Hvl^ ^ILl(0TT°°(f^)) 

— ^lk<^llL°°(0,T;L2(Q;R™xn)) ||c(^£)\/k£| + l^£l||Too(o,T;L2(C2;Rnxn))||'\/lC£ ~ ^£l|lLl(0,r;L°°(Q)) 

— ^lk^llL°°(0,T;L2(Q;Rnxn)) ||k(^i£)| 5(2:£)||j;^cx>(o,T;Ll(f^;*”^")) “ ^sllLl(0TT°°(f2)) 

— ^lk<^IL°°(0,T;L2(n;RriXn)) (^||bD {Ce, e{Ue), Ze)\\ + |i C£ || (o,T;I,2(Q;Rnxn))) X 


X 


11^^ 


'£NlLl(0,r;L°°(Q))' 


(23) 
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By using the a priori estimates in Lemma 13.21 (i) and (v), we obtain 


ll'^s|L°°(0,T;L2(r2;R'iX'*)) (ll^ IIL°°(o^T;LlIIIIL°°(0,r;L2(f7;]Rnxn))) < C 

and by using the convergence properties (HEl) and (IM and the compact embedding 
C{VIt), we end up with 

llVlCe — ^£||Li(0,T;L°°(C2)) ^ 

as e —>■ 0"^. Therefore, the right hand side of ()23p converges to 0 as e —)■ O"*". 

• Exploiting the boundedness properties in Lemma 13.21 (i), (iii) and (iv) and the convergence 
properties (I19cp and (I21ap . the third integral on the right hand side of (1221) converges to 0 as 
e ^ 0+. 


We have shown limsupj,_^Q+= 0. Futhermore, the boundedness of {V^g} in LP(n'r;M”’) (see 
Lemma [32]) and the convergence property (I 21 ap show Bg —>• 0 as e —> O"''. 

We conclude with limj_^o+ — Vzl^d^dt = 0. Therefore, (b) is also shown. 

Property (c) is an immediate consequence of the dehnition of and the convergence property 
6 ° ^ 6 ° in 

To prove (d), we define := {z^ > 0} D {z > 0}. Consequently, we get 






(24) 


and the convergence 


Iwe —>■ l{^>o} pointwise in Qt and in for all <7 > 1 (25) 

as e —>■ 0 + by using ^ z in Ot- Calculating the weak L^( 07 ’; M”^”')-limits in (IMl) for e —>■ O"'' 
on both sides by using the already proven convergence properties, we obtain W^eic,e, z)l^z>o} = 
tCel{ 2 >o}- The remaining convergence properties in (d) follow from Lemma 13.21 fvl. □ 


Corollary 3.4 Let t G [0,T] and U CC {z{t) = 0} be an open subset. Then U C {^^(s) > 0} for 
all s G [0, t] provided that e > 0 is sufficiently small. More precisely, there exist 0 < rj such that 

■Ze(s) >rjinU 


for all s G [0, t] and for all 0 < e ^ 1. 

Proof. We have U x [0, T] C {z > 0} and infj-^ i)gI7x[o t] s) > 0 by 2 : G C{LIt)- The claim follows 
from Zc ^ z in C{LIt) according to Lemma 13.31 (b). □ 

To obtain a-priori estimates for the chemical potentials {ne} in a local sense, we make use of 
the so-called conical Poincare inequality for star-shaped domains cited below. 


14 










Theorem 3.5 (Conical Poincare inequality jBK98]) Suppose that Vl <ZW^ is a hounded and 
star-shaped domain, r > 0 and 1 < p < oo. Then there exists a constant C = C{Q,p,r) > 0 such 
that 

f |i(;(x) — icq dx < C f |Vic(x)|^(5^(x) dx 

Jn ’ Jn 

for all w G (7^(0), where the 5'^-weight is given by 

'■= / w{x)6^{x)dx, 5(x) := dist(x, 9n). 

Jn 

By a density argument, the statement is, of course, also true for all w G 

Lemma 3.6 (A-priori estimates for 

(i) Interior estimate. For every t G [0, T] and for every open cube Q CC {z{t) > 0} H hi, there 


exists a C > 0 such that for a// 0 < e <C 1 

ll/^e||L2(o,f;_ffi(Q)) < C. 


(26) 


(a) Estimate at the boundary. For every t G [0, T] and every xq G {z{t) > 0} H dQ, there exist a 
neighborhood U of xq and a C > 0 such that for a// 0 < e -C 1 


\\Ts\\L-^(Q,t-,m{unn)) < C. 


(27) 


Proof. 


(i) Let t G [0, T] and Q CC {z{t) > 0} n 0 be an open cube. We consider the Lipschitz domain 
Q := Brj{Q) := {x G M" |dist(x,(5) < e}, where 77 > 0 is chosen so small such that Q CC 
{z{t) > 0} n n. We define the following function 


C(x) 


A(x) := dist(x, dQ) if x G Q, 
0 else. 


is a Lipschitz function on hi with Lipschitz constant 1. By Rademacher’s theorem (see, for 
instance, |Zie89[ Theorem 2.2.1]), f is in W^’°°(il). Now, we can test (fT^ with (. Then, by 
using the previously proven a-priori estimates, we obtain boundedness of 

Psix, s)A(x) dx < C (28) 

Jq 

with respect to a.e. s G (0, T) and e. 

Furthermore, there exists an r/ > 0 such that zds) > rj in Q for all s G [0,t] and for all 
0 < e ^ 1 (see Corollary I3.4li . Thus, by assumption ([3]), m^izs{s)) > t]' > 0 holds in Q for 
all s G [0, t] and all 0 < e <C 1 for a common constant r]' > 0. Consequently, we get by the 
a-priori estimate for ps 

ll^/^£llL2(Qx[0,t]) - ^ 

for all e. Applying Theorem 13.51 (we put inFt, = Q,r=p=2 and w = Peis) for s G [0, tj), 
integrating from 0 to t in time and using boundedness properties (I25|) and (12^ . we obtain 
boundedness of \\P£M\L^(Qx[ot]) boundedness of \\p£\\L^{Qx[o,t]) with respect to 0 < 

e -C 1. Together with (1291) . we get the claim (1^ . 
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(ii) Since {z{t) = 0} C 0 is a closed set, we can find a neighborhood U CW^ \ {z{t) = 0} of xq. 
Furthermore, since hi has a C^-boundary, there exists a C^-diffeomorphism vr : (—I,!)” —>■ U 
with the properties 

. 7r((-l,ir-i X (-1,0)) C5d, 

• 7r((-i,i)”-i X { 0 }) c on, 

• 7r((-l, 1)^-1 X (0,1)) C \ a 

Let : (—1,1)" —>■ (—1,1)"' denote the reflection x i->- (xi,..., x„_i, —Xn) and T := 7ro-do7r“^. 
Furthermore, let /ig G L‘^{0,t] (U)) be dehned by 


fie{x,s) 


lle{x, s) if X € U n fi, 

^s{T{x),s) if X G t/\ id. 


Let Q CC [/ be a non-empty open cube with xq G Q. Then, integration by substitution with 
respect to the transformation T yields for a.e. s G (0, t) 

/ ^^(x, s)A(x) dx = / /i£(x, s)A(x) dx 

Jq Jqnn 

+ [ /Xe(x, s)A(T(x))| det(VT(x))| dx, (30) 

where the Lipschitz function A : —)• M is given by 


A(x) 


dist(x, dQ) if X G Q, 

0 if X G \ Q. 


We are going to show that both terms on the right hand side of (f3(l are bounded with respect 
to e and a.e. s G (0,t). 

• Testing (fl^ with the function C = -^ yields 




(Vc,(s)-VA + ^,,(c,(s))A) dx 


+ [ W^,c(Ce('S)T(w(s)),2;e(s))Adx. 

jQnn 

By the already know a-priori estimates, every integral term on the right hand side is 
bounded w.r.t. £ and a.e. s G (0, T). 

• The function 


C(x) 


(A(r(x)))| det(vr(x))| if X G T{Q \ O), 

0 if X G \ T{Q \ n) 


is a Lipschitz function in because: 


(31) 
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- A o T is a Lipschitz function in [/ D and A o T = 0 in n O) \ T(Q \ ri). The first 
property follows from the Lipschitz continuity of A and of T (note that T is a 
diffeomorphism). The latter property can be seen as follows. Assume the contrary. 
Then, we find an x € (17 fl fl) \ T((5 \ fl) such that A(T(x)) > 0. By the definition 
of A, we get T(x) € Q. Since x G hi, it follows T{x) 0 hi by the construction 
of T. Therefore, T(x) G Q \il. This gives x = T(T(x)) G T(Q \ n) which is a 
contradiction. 

- I det(VT)| is a Lipschitz function in [/ n n (note that T is a C^-diffeomorphism). 

Testing (fl^ with ( from (l3Tl) yields 


lT{Q\n) 


/x,(s)(Aor)|det(Vr)|dx 


lT{Q\n) 


Vce{s) • V((A o T)\ det(Vr)l) dx 


+ / (T,c(ce(s))+ lT,c(c£(s),e(rie(s)),Ze(s)))(AoT)|det(VT)|dx. 

JTiQ\n) 

By the already know a-priori estimates, every integral term on the right hand side is 
bounded w.r.t. e and a.e. s G (0, t). 

For VJlsis), we also get by integration via substitution: 

[ [ |V//£(x,s)pdxds 
Jo Jq 

[ [ |V/i£(x,s)pdxds+ [ [ |V/re(T(x),s)p|VT(x)pdxds 

Jo JQnn 

[ [ |V/r£(x,s)pdxds 

Jo Jonn 


< 


0 JQ\n 


+ f [ |v/i,(x,s)nvr(r(x))ndet(vr(x))|dxds. 

Jo JT(Q\n) 


(32) 


Since Q fl fl CC {z{t) > 0} and T{Q \ CC {z{t) > 0}, we deduce Zs{s) > rj on Q CiQ and 
on T{Q\ ri) for all s G [0,t] and for all sufficiently small 0 < e (see Corollary [3lT|) . Thus, 
VjUs is bounded in L^{{Q fl O) x (0, t);]R"') and in Lp‘{T{Q \ Q) x (0,t);M") with respect to 
0 < e <C 1 by also using the a-priori estimate for the property m G C([0,1]; M+) 

and assumption ([3]). 

Therefore, the left hand side of (I32p is also bounded for all 0 < e <C 1. The Conical Poincare 
inequality in Theorem 13.51 yields boundedness of ^^A in L?‘{Q x (0,t)). Finally, we can find a 
neighborhood B C Q of xq such that //g is bounded in L‘^{0,t; H^{V)). □ 


Due to the a-priori estimates for {ne} and {u^}, the limit functions and u can only be expected to 
be in some space-time local Sobolev space (see Section [2]). In the sequel, it will be necessary 

to represent the maximal admissible subset of the not completely damaged area, i.e. ^£){{z > 0}), 
as a union of Lipschitz domains which are connected to parts of the Dirichlet boundary D. To 
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this end, we define by F the shrinking set F := {z > 0} until the end of this work and obtain the 
following result. 

Lemma 3.7 (cf. [HK12L Lemma 4.18]) There exists a function u G M”) such 

that e(tt) = e a.e. in ^r{F) and u = b on the boundary Dt n 2tr(L"). 

A related result can be shown for the sequence {ne} by exploiting the estimates in Lemma 13.61 To 
proceed, we recall a definition introduced in |HK12] . 

Definition 3.8 (cf. |HK12L Definition 4.1]) Let H F Ll be a relatively open subset. We call a 
countable family {Uk] of open sets Uk CC H a fine representation for H if for every x & H there 
exist an open set U C with x € U and an k gN such that U HQ C Uk- 

Lemma 3.9 Let a sequence {tm} Q [0, T] containing T be dense. There exist a fine representation 
{UJfi }k&N for F{tm) for every m gN, a function n G L^H^ ^oci^) ® subsequence of {fie} (also 
denoted by {fie}) such that for all k,m Gf^ 

fie^ fi mL\0,tm;H\U]fi)) (33) 


as e —)• O’*". 

Proof. A fine representation {C/™}fceN of F{tm) can be constructed by countably many open cubes 
Q CC F{tjn) n n together with finitely many open sets of the form U DLl such that Lf satisfies (|27p 
from Lemma 13.61 (ii). Then, we have for each k,m gN the estimate 

for all 0 < £ -C 1. By successively choosing sub-sequences and by a diagonal argument, we obtain 
a fi G LIH). iq(,(T’) such that (l33]) is satisfied (cf. proof of |HK12l Lemma 4.18]). □ 

The a-priori estimates and the convergence properties of {ze} in Lemma 13.31 and of {fie} in 
Lemma 13.91 respectively, yield the following corollary. 

Corollary 3.10 It holds for £ —)■ O'*'." 

m{ze)Vfie m{z)Vfi in L^(F;M"), 

m{ze)Vfie 0 in Lfi{LlT \ F;]R"'). 

Now, we have all necessary convergence properties to perform the degenerate limit in (|inp - (|15p . To 
proceed, we need the following auxiliary result. 

Lemma 3.11 Let {t^} and {Ufi^} be as in Lemma I,?..91 Then, for every compact subset K F F 
there exist a finite set I Ffl, values mk GN, k G I, and functions ipk S C°°{LIt), k G I, such that 

(i) 

(a) supp('0fc) c L'™'' X 
Dkei'^k = 1 on K. 
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Proof. We extend the family of open sets {14™} given by I4™ := UJP x (0, in the following way. 
De fin e 

T’ ■= {{Wnk,mm I WJP C is open with Wf n Dt = x (0,t™j}. 

We see that V is non-empty and that every totally ordered subset of V has an upper bound with 
respect to the ” < “ ordering defined by 

{WfP} < {W^} ^ C W™ for all A:, m G N. 

By Zorn’s lemma, we find a maximal element {14”}- It holds 


rc U 

/c,m,EN 


m 

k * 


(34) 


Assume that this condition fails. Then, because of F n Qt = Ufc meN there exists a p = (x, t) G 

F n 9(Dr) with p 0 (Jfc.meN _ 

Let us consider the case t < T. Since F C Dy is relatively open, we find an niQ G N with 
X G F{tmo) and > L By the fine representation property of {A/™°}fceN for F(tmo), 
open set U C M"’ with x € U and A:o G N such that [/ n D C U^°. 

The family {W™} given by 


WF:= 


14™ U t/ X (—00, tmo) if k = ko and m = mo, 

else, 


for 


satisfies {Wf} G F and p G Ufc,meN W™ which contradicts the maximality property of {Vjf^}. 

In the case t = T, we also hnd ko,mo G N and an open set U C M” with x G 17 such that 
17 n D C and tmo = T. The family {W™} given by 


WF ■■= 


14™ U 17 X 

VF 


if k = ko and m = mo, 
else, 


also contradicts the maximality of {14™}- Therefore, ()34p is proven. 

Heine-Borel theorem yields 

Kc[J v;'" 

k&I 

for a finite set / C N and values G N, A: G /. Together with a partition of unity argument, we 
get functions ifk € C°°{Qt) such that (i)-(iii) hold. □ 

The degenerate limit e ^ O'*' can be performed as follows: 

• We define the strain by e := ejj? G and obtain for the remaining variables 

c€ L^{0,T-,H\n))nH\0,T-,iH\n)r), uGLjHlU^DiFfR^), 
z G L°-{0,T;W^’P{n))nH\0,T-,L^{n)), p G 


with e = e(tt) in 2t(F). 
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• Passing to the limit e —> O’*" in (fTOjl , (fT3|) and (fTHIl imply properties ([5]) , ([8]) and ([9]) , cf. |HK12| . 

• Usine Lemma l,S.31 lal and Corollary 13.101 we can oass to —>■ O'*' in (ffH) and obtain Q. 

• Let C £ -b^(0, T; with supp(C) C F be a test-function. Furthermore, let {^pl} be a 

partition of unity of the compact set K := supp(C) according to Lemma 13.111 For each I £ N, 
we obtain snpp{C'4’i) ^ x [0, tmj- Then, integrating (fT^ in time from 0 to tmi, testing 
the result with (ipi and passing to e —>• 0'*' by using Lemma [T3l and Lemma iTOl show 



dx ds 



Vc • V(CV'z) + ^,c(c)C^i + bF,c(c, e, z)C'<pi'^ 


dx ds. 


Summing with respect to I G I and noticing = 1 on supp(C) yield ([7]). 


In conclusion, the limit procedure in this section yields functions (c, e, u, z, /i) with e = e(u) in 
2l£)(F) which satisfy properties (ii)-(v) of Definition 12.11 In particular, the damage function 2 has 
no jumps with respect to time. However, we cannot ensure that {z > 0} equals 21 d({^ > 0}) and, 
moreover, if F \ 21/)({z > 0}) 7^ 0, it is not clear whether u can be extended to a function on F 
such that e = e(u) also holds in F. This issue is addressed in the next section where such limit 
functions are concatenated in order to obtain global-in-time weak solutions with hneness t/ > 0 by 
Zorn’s lemma. 


4 Existence results 

In this section, we are going to prove the main results of this paper. 

Theorem 4.1 (Maximal local-in-time ^veak solutions) 

Let b G c° G H^{n) and z° eW^’P(n) with 0 < k < z° < 1 in Q be 

initial-boundary data. Then there exist a maximal value T > 0 with T < T and functions c, u, z, 

/i defined on the time interval [0,T] such that {c,u,z,n) is a weak solution according to Definition 
\2.1[ Therefore, ifT <T, {c,u,z,fi) cannot be extended to a weak solution on [0,T -|-e]. 

Proof. Zorn’s lemma can be applied to the set 

F := {(T, c,u, z, fj,) \ 0 < T <T and (c, u, z, /i) is a weak solution on 
[0, r] according to Definition 12.ll) 

to hnd a maximal element with respect to the following partial ordering 

(Ti, c\, Ur, z \, ) F (T2 , C2, U2 , Z2 ,1^2) Tr F T2 , C 21 j ? ^2 1 [^0 Ti] ? 

^2|[o,Ti] ~ = hi- (35) 

Indeed, F / 0 by the result in Section[3J More precisely, since z G L°°(0, F; lT^’P(D))nF^(0, T; L^(D)) 
and since 0 < k < 2;*^, we find an e > 0 such that {z{t) > 0} = ^jo{{z{t) > 0}) for all t G [0,e]. For 
the proof that every totally ordered subset of F has an upper bound, we refer to |HK12j . □ 
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The proof of global-in-time existence of weak solutions with fineness r] requires a concatenation 
property (see Lemma I4.4p which is, in turn, based on some deeper insights into the T-limit (£ 
introduced in Section [3j To this end, it is necessary to have more information about the recovery 

p 

sequences for t d- 

We will introduce the following substitution method: Assume that u G minimizes 

e(-), z) with Dirichlet data ^ on D. Then, by expressing the elastic energy density W in terms 
of its derivative he. 

W(c,e,z) = ^We(c, e, z) : e + ^zip^(c) : e + z(p^(c}, 

and by testing the momentum balance equation with Q = u — u for a function u G M”) with 

u = ^ on D, the elastic energy term in Tg can be rewritten as 

[ We{c,e{u),z)dx = [ {g{z) + e) ((p^e{u) : e{u)+ ]-(p‘^{c) : {e{u) + e{u)) + dx. (36) 

Jn Jn \ z 

For convenience, in the following proof, we define the density function W as 

if^e : ei -h : (e -h ei) -t- 



We{c,e,ei,z) := {g{z) -he) 


Lemma 4.2 For every c G ^ G VF^’°°(H) and z G kF^’^(H) there exists a sequence 5g O'*' 

p 

such that (c, {z — hg)"'') —>■ {c,^,z) is a recovery sequence for —> 5- 

Proof. We follow the idea of the proof in |HK12[ Lemma 4.9], But here we have to deal with the 
additional concentration variable which complicates the calculation. Let {cg,^g,Zg) —>■ (c,^,z) be a 
recovery sequence and 6g —)■ 0+ such that {z — 6g)~^ < Zg. Consider 


deic,^, {z - (5^) + ) -ds{Ce,^s,Ze) = deic,f,, (z - 4) + ) - deic, f,, Zg) + ^g{c, Zg) - de{Ce,^e,Ze) 


Be 


We have A,. < 0 since Fg{c,e{f, -|-C)) (-2 — (^e)^) < -l-C)!-^;^) for all C, G Now, we 

focus on the second term of the right hand side. Let Ug,Vg G Lf|j(0;M”) be given by 


Ug= argmin ^^(c, e(^-I-C), ^e), Vg = argmin Fg{cg,e{^g + (), Zg). 

By using the substitution (l36|) for Wg{c.,e{f, + Ug),Zg) with test-function u = Vg and (l36l) for 
Wg{cg, e(^e -h Vg), Zg) with test-function u = Ug, we obtain a calculation as follows: 


:(c, c(C T rtg), 2^) J'zipg, e{fg Vg), Zg,) 

= (f^(c,e(? + w),e(? + Ve),Zs + £) -W{cg,e{fg + Vg),e{fg + Ug),Zg + £)'^ dx 

= Jj^size) + e) -h Ug) : e(^ -h Vg) - (p^€{f,g + Ug) : €{f,g + Vg) 

+ ■ ^(2^ + Ue + Vg) - ]^^p^{cg) : e(2^ + Vg + Ug) -h ip^{c) - (f^icg)'^ dx 
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= jj-Size) + e) : e(0 - ■ 4Ce) + + Ve) : e(C - Ce) 

+ <y?^(c) : e(C - Ce) + - ^^(Ce)) ■ €{2^e + Ue + Ve) + (f^ic) - (f^iCe)^ dx 

< Jj^size) + e) : e(0 “ V^^e(6) : e(6) + ^^(c) ■ e(? - ?e) + ^^{c) - dx 

+ \\{g{Ze) + e)(p^e{Ue + ^^£)||L2(t^)||e(^ - 6 )||l2(q) 

+ 2 “ ¥^^(Ce)llL 2 (r 2 ) (11(5(2^6) + ^)e(?£ + ^£)lll,2(n) + \\{g{Ze) + e)e(Ce + 't’£:)llL 2 ( 0 )) 

Using the convergence properties ^ c in H^{Q), ^ ^ in Ty^’°°(r2), 2 :^ ^ z in W^’P{Q) 

and the boundedness of J>(c, e(^ + Ue),Ze) and Te{ce,e{^e + Ve)-,Ze) with respect to e, we conclude 
lim supj,_j.Q+-Be < 0. Now we can proceed as in the proof of |HK12[ Lemma 4.9] and the claim 
follows. □ 


r 

Remark 4.3 The knowledge of such recovery sequences for d gives also more information 

p p 

about (B. In particular, we obtain that the result in Lemma If-^ with <Be <B instead of ^e d 
also holds true and, moreover, that the following properties are satisfied (cf. IHK1S\. Corollary 4-10, 
Lemma 4-ilk])- 

• (B{c,c,1fz) <<Bic,tz) ycGH^{n),vc€W^^^{n-w^),yz€W^’P{n) 

VB C n open with Ipz € W^'P{Ll), 

• B{c,i,z) <£{c,e{u),z) ^c ^ ^z with 

0 < z < 1, Vtt G HI^^{{z > 0}; M”) with u = f, on D Ci {z > 0}. 

Lemma 4.4 Let ti < t 2 < ts be real numbers and let p > 0. Suppose that 
q := {c,e,u,z,Jl,F) is a weak solution with fineness p on [ti,t2]) 

q := {c,e,u,z,fi,F) is a weak solution with fineness p on [t 2 ,t 3 ] withTf^ = C:(c(t2), ^(^2);^^(^2)) 

(the value for the weak solution q in Definition \2. 1\) . 

Furthermore, suppose the compatibility condition 0(^2) = c(t2) and z^(t 2 ) = z~it2)^2ijj({z-{t 2 )>o}) 
and the Dirichlet boundary data b G 

Then, we obtain that q := {c,e,u,z,fi,F) defined as q\[t-i,t 2 ) ■— Q g\[t 2 ,t 3 ] ■= q is a weak 
solution with fineness p on [tijts]. 

Proof. Because of the properties in Remark 14.31 we can prove the following crucial energy estimate 
at time point t 2 '. 

lim essinf £’(c(r), e(r), 2;(r)) = lim ess inf £’(c(r), e(r), (r)) 

s^t~ re{s,t2) s^t~ Te{s,t2) 

> lim_ essinf £'(c(r),e(u(r)),2;“(r)l2i3({^-(^)>0})) 

S—>-t 2 '^^(^ 5 ^ 2 ) 

> lim_ essinf (£(c(r),6(T),2;“(r)l2i^({^-(^)>0})) 

S—>-t 2 '^€(^ 5 ^ 2 ) 
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> (£(c(i2),6(t2),x) 

> >^{c{t2),b{t2),z'^{t2)) 


with X := -2 (i2)lf| j )2lD({^“(r)>o})- With this estimate, we can verify the claim by the same 
argumentation as for |HK12[ Lemma 4.21], □ 


Theorem 4.5 (Global-in-time weak solutions with fineness rj) 

Let b G c° € H^{n) and G W^’P{n) with 0 < < 1 in Cl and {2° > 0} 

admissible with respect to D be initial-boundary data. Furthermore, let rj > 0. Then there exists a 
weak solution {c,e,u,z,n) with fineness i] > 0 according to Definition \2.S\. 

Proof. This result can also be proven by using Zorn’s lemma on the set 

V := {{f,c,e ,u, z, fj., F) \ 0 < T < T and (c, e, u, z, g,, F) is a weak solution with fineness rj on 

[0,T] according to Definition 12.31} 

with an ordering analogously to (1351) . The assumptions for Zorn’s lemma can be proven as in The¬ 
orem 10 (see |HK12[ Proof of Theorem 4.1]). To show that a maximal element from V is actually 
a weak solution with fineness p on the time-interval [0, T], we need the concatenation property in 
Lemma 01 Indeed, if a maximal element q is only defined on a time-interval [0,T] with T < T 
we can apply the degenerated limit procedure in Section [3] to the initial values c(T) and z{T) to 
obtain a new limit function q. By exploiting Lemma [4.41 5 is a weak solution with fineness r] on the 
time-interval [0, T -\- e] for a small e > 0 which contradicts the maximality of q. □ 


References 


[BB99] J.W. Barrett and J.F. Blowey. Finite element approximation of the Cahn-Hilliard equation with concen¬ 
tration dependent mobility. Math, of Computation, 68(226):487-517, 1999. 

[BCD"*'02] E. Bonetti, P. Colli, W. Dreyer, G. Gilardi, G. Schimperna, and J. Sprekels. On a model for phase 
separation in binary alloys driven by mechanical effects. Physica D, 165:48-65, 2002. 

[BK98] S. M. Buckley and P. Koskela. New Poincare inequalities from old. Annales Academice Scientiarum 
FenniccE Mathematica, 23:251-260, 1998. 

[BMIO] S. Bartels and R. Muller. A posteriori error controlled local resolution of evolving interfaces for generalized 
Cahn-Hilliard equations. Interfaces and Free Boundaries, 12(l):45-73, 2010. 

[BMR09] G. Bouchitte, A. Mielke, and T. Roubfcek. A complete-damage problem at small strains. ZAMP Z. Angew. 
Math. Phys., 60:205-236, 2009. 

[BP05] L. Bartkowiak and I. Pawlow. The Cahn-Hilliard-Gurtin system coupled with elasticity. Control and 
Cybernetics, 34:1005-1043, 2005. 

[BS04] E. Bonetti and G. Schimperna. Local existence for Fremond’s model of damage in elastic materials. Contin. 
Mech. Thermodyn., 16(4):319-335, 2004. 

[GH58] J.W. Cahn and J.E. Hiliard. Free energy of a uniform system, i. interfacial free energy. Journal of Chemical 
Physics, 28:258-267, 1958. 

[GMPOO] M. Garrive, A. Miranville, and A. Pietrus. The Gahn-Hilliard equation for deformable elastic continua. 
Adv. Math. Set. AppL, 10(2):539-569, 2000. 

[DM01] W. Dreyer and W.H. Muller. Modeling diffusional coarsening in eutectic tinlead solders: A quantitative 
approach. Intemat. J. Solids and Structures, 38(8): 1433-1458, 2001. 


23 






[EG96] 

[E1189] 

[FK09] 

[GarOO] 

[GarOSa] 

[Gar05b] 

[Gia05] 

[GRWOl] 

[GUE+07] 

[Gur89] 

[HCW91] 

[HKll] 

[HK12] 

[HK13] 

[KRZll] 

[Mer05] 

[Miell] 

[MR06] 

[MRZIO] 

[MTIO] 

[PZ08] 

[RR12] 

[Sim86] 

[USG07] 

[WeiOl] 

[Zie89] 


C.M. Elliott and H. Garcke. On the Cahn-Hilliard equation with degenerate mobility. SIAM J. Math. 
Anal, 27:404-423, 1996. 

C.M. Elliott. The Cahn-Hilliard Model for the Kinetics of Phase Separation. In Mathematical Models for 
Phase Change Problems (ed. J.-F. Rodrigues), volume 88 of International Series of Numerical Mathematics, 
pages 35-73. Birkhauser Basel, 1989. 

J. R. Fernandez and K. L. Kuttler. An existence and uniqueness result for an elasto-piezoelectric problem 
with damage. Math. Mod. Meth. Appl. Sci., 19(l):31-50, 2009. 

H. Garcke. On mathematical models for phase separation in elastically stressed solids. Habilitation thesis. 
University Bonn, 2000. 

H. Garcke. Mechanical Effects in the Cahn-Hilliard Model: A Review on Mathematical Results. In 
A. Miranville, editor. Mathematical Methods and Models in phase transitions, pages 43-77. Nova Science 
Publ., 2005. 

H. Garcke. On a Cahn-Hilliard model for phase separation with elastic misfit. Annales de I’Institut Henri 
Poincare (C) Non Linear Analysis, 22(2):165 - 185, 2005. 

A. Giacomini. Ambrosio-Tortorelli approximation of quasi-static evolution of brittle fractures. Calc. Var. 
Partial Differ. Equ., 22(2):129-172, 2005. 

H. Garcke, M. Rumpf, and U. Weikard. The Cahn-Hilliard equation with elasticity: Finite element 
approximation and qualitative studies. Interfaces Free Bound., 3:101-118, 2001. 

M.G.D. Geers, R.L.J.M. Ubachs, M. Erinc, M.A. Matin, P.J.G. Schreurs, and W.P. Vellinga. Multiscale 
Analysis of Microstructura Evolution and Degradation in Solder Alloys. Internatilnal Journal for Multiscale 
Computational Engineering, 5(2):93-103, 2007. 

M.E. Gurtin. On a Nonequilibrium Thermodynamics of Capillarity and Phase. Quarterly of Applied 
Mathematics, 47(1):129-145, 1989. 

P. G. Harris, K. S Chaggar, and M. A. Whitmore. The Effect of Ageing on the Microstructure of 60:40 
Tin-lead Solders. Soldering & Surface Mount Technology Improved physical understanding of intermittent 
failure in continuous, 3:20-33, 1991. 

C. Heinemann and C. Kraus. Existence of weak solutions for Cahn-Hilliard systems coupled with elasticity 
and damage. Adv. Math. Sci. Appl., 21(2):321-359, 2011. 

C. Heinemann and C. Kraus. Complete damage in linear elastic materials — modeling, weak formulation 
and existence results. WIAS preprint no. 1722 (submitted), 2012. 

C. Heinemann and C. Kraus. Existence results for diffuse interface models describing phase separation 
and damage. Eur. J. Appl. Math., 24(2):179-211, 2013. 

D. Knees, R. Rossi, and C. Zanini. A vanishing viscosity approach to a rate-independent damage model. 
WIAS preprint no. 1633. WIAS, 2011. 

T. Merkle. The Cahn-Larche system: A model for spinodal decomposition in eutectic solder; modelling, 
analysis and simulation. PhD-thesis, Universitat Stuttgart, Stuttgart, 2005. 

A. Mielke. Complete-damage evolution based on energies and stresses. Discrete Contin. Dyn. Syst., Ser. 
S, 4(2):423-439, 2011. 

A. Mielke and T. Roubfcek. Rate-independent damage processes in nonlinear elasticity. Mathematical 
Models and Methods in Applied Sciences, 16:177-209, 2006. 

A. Mielke, T. Roubfcek, and J. Zeman. Complete Damage in elastic and viscoelastic media. Comput. 
Methods Appl. Mech. Engrg, 199:1242-1253, 2010. 

A. Mielke and M. Thomas. Damage of nonlinearly elastic materials at small strain — Existence and 
regularity results. ZAMM Z. Angew. Math. Mech, 90:88-112, 2010. 

I. Pawlow and W. M. Zaj§,czkowski. Measure-valued solutions of a heterogeneous Cahn-Hilliard system in 
elastic solids. Colloquium Mathematicum, 112 No.2, 2008. 

E. Rocca and R. Rossi. A degenerating PDE system for phase transitions and damage. arXiv:1205.3578vl, 

2012 . 

J. Simon. Compact sets in the space {0,T-, B). Annali di Matematica Pura ed Applicata, 146:65-96, 
1986. 

R.L.J.J. Ubachs, P.J.G. Schreurs, and M.G.D. Geers. Elasto-viscoplastic nonlocal damage modelling of 
thermal fatigue in anisotropic lead-free solder. Mechanics of Materials, 39:685-701, 2007. 

U. Weikard. Numerische Losungen der Cahn-Hilliard-Gleichung und der Cahn-LarchAGleichung. PhD- 
thesis, Universitat Bonn, Bonn, 2001. 

W.P. Ziemer. Weakly differentiable functions. Springer-Verlag New York, Inc., 1989. 


24 


